Supersymmetry breaking as a quantum phase transition 
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We explore supersymmetry breaking in the light of a rich fixed-point structure of two-dimensional 
supersymmetric Wess-Zumino models with one supercharge using the functional renormalization 
group (RG) . We relate the dynamical breaking of supersymmetry to an RG relevant control param- 
eter of the superpotential which is a common relevant direction of all fixed points of the system. 
Supersymmetry breaking can thus be understood as a quantum phase transition analogous to similar 
transitions in correlated fermion systems. Supersymmetry gives rise to a new superscaling relation 
between the critical exponent associated with the control parameter and the anomalous dimension 
of the field - a scaling relation which is not known in standard spin systems. 
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Supersymmetry is an important guiding principle for 
the construction of models beyond the standard model 
of particle physics, as it helps parametrizing the problem 
of a large hierarchy of scales, supports the unification 
of gauge couplings, and facilitates an intricate combina- 
tion of internal symmetries with the Poincare group. On 
the other hand, supersymmetry must be broken in na- 
ture in a specific way without violating these attractive 
features and other constraints imposed by the standard- 
model precision data. In addition to a variety of break- 
ing mechanisms often involving further hidden gauge and 
particle sectors, dynamical supersymmmetry breaking py 
can play a decisive role in the formation of the low-energy 
standard model. Generically, symmetry breaking is re- 
lated to collective condensation phenomena, often requir- 
ing a nonperturbative understanding of the fluctuation- 
induced dynamics. 

From the viewpoint of statistical physics, symmetry 
breaking and phase transitions between ground states 
of a differently realized symmetry are often related to 
fixed points of the renormalization group (RG). In ad- 
dition, the quantitative properties of the system near a 
phase transition are determined by the RG flow in the 
fixed-point regime. As the understanding of symmetry 
breaking phenomena in quantum field theory and parti- 
cle physics has often profited from analogies to statistical 
systems, this paper explores supersymmetry breaking in 
the light of phase transitions and critical phenomena. 

As a simple example, we concentrate on Wess-Zumino 
models with Af — 1 supercharges in two dimensions in 
Euclidean spacetime. Dynamical supersymmetry break- 
ing in this and many other supersymmetric systems is 
governed by a control parameter of the superpotential 
which plays the role of a bosonic mass term. This is very 
reminiscent to phase transitions in strongly correlated 
fermionic systems: here, symmetry breaking is governed 
by a composite bosonic order parameter, e.g., a super- 
fluid condensate, the mass term of which serves as the 
control parameter of the phase transition. For instance, 
the inverse interaction strength of attractively interact- 
ing electrons on a honeycomb lattice governs the phase 



transition between the semimetallic and the superfluid 
phase as discussed in the context of ultracold atoms 
or graphene ^] . As this transition occurs at zero temper- 
ature, the critical value of the control parameter marks 
a quantum critical point of a quantum phase transition. 

A standard approach to quantum phase transitions in 
strongly correlated fermion systems is the Hertz-Millis 
theory [1, \^ where composite bosonic fields are intro- 
duced by a Hubbard- Stratonovich transformation. Sub- 
sequently, the fermionic degrees of freedom are integrated 
out, leaving a purely bosonic description which is dealt 
with in a local approximation. As this strategy is not 
always sufficient, a treatment of fundamental fermionic 
and composite bosonic degrees of freedom on equal foot- 
ing beyond the Hertz-Millis theory has recently proved 
successful @. 

Now, supersymmetric systems by definition are com- 
bined bosonic-fermionic systems with a high degree of 
symmetry. A treatment of these degrees of freedom on 
equal footing is even mandatory to preserve supersym- 
metry manifestly. Apart from this symmetry, we will 
show that supersymmetry breaking has many similari- 
ties to a quantum phase transition in strongly correlated 
fermion systems. Beyond these similarities, supersymme- 
try invokes additional structures which go beyond those 
known in statistical physics. 

Euclidean two-dimensional J\f = 1 Wess-Zumino mod- 
els can be defined by the Lagrange density in the off-shell 
formulation C = £kin + 'Cint with 

/:ki„ ^^d^cj.d^'c^ + - '-d.i^r^ ~ If', (i) 



(2) 



with superpotential W {(()), bosonic field (f), fermionic 
fields V' and V' and auxiliary field F. The prime denotes 
the derivative with respect to 0, and we use a chiral rep- 
resentation for the 7 matrices. This model allows for 
dynamical supersymmetry breaking if the highest power 
of the bare superpotential is odd. 

We use the functional RG to calculate the effective 
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action T of the theory. The functional RG can be formu- 
lated as a flow equation for the effective average action Tk 
• This is a scale-dependent action functional which in- 
terpolates between the classical action S = J (^kin+'Cint) 
and the full quantum effective action F. The interpola- 
tion scale k denotes an infrared (IR) regulator scale, such 
that no fluctuations are included for fc — > A (with A be- 
ing the UV cutoff), implying F^^a — > S. For fc ^ 0, all 
fluctuations are taken into account such that Tk^o F 
is the full quantum solution. The effective average action 
can be determined from the Wetterich equation [7j] 



dtVk = \ STr 



.(2) 



dfRk 



t = In- 



A' 



(3) 



which defines an RG flow trajectory in theory space 
with S serving as initial condition. Here, (F^'^M = 

V / ab 

j^Ffcj^ where the indices a,b summarize field com- 
ponents, internal and Lorentz indices, as well as space- 
time or momentum coordinates, i.e., 4'"'" — {4),F,'ip,^). 
The regulator function Rk is introduced in the form of 
a quadratic contribution to the classical action, AS'fe — 
J'^^Rk'^. Rk can thus be viewed as a momentum- 
dependent mass term, implementing the IR suppression 
of modes below k. Different functional forms of Rk cor- 
respond to different RG schemes. Physical quantities do 
not depend on the regulator. 

For a manifestly supersymmetric flow equation, we de- 
rive the regulator from a D term of a quadratic superfield 
operator, similarly to our construction for supersymmet- 
ric quantum mechanics [s^. The most general supersym- 
metric cutoff action in the off-shell formulation in mo- 
mentum space is given by 



ASk~- J d p(rii/)7*-0 — 2ri0F-fp r20 +r2'>pf4> — r2F-). 

The regulator shape functions ri^2 = ?'i,2(p^) determine 
the precise form of the Wilsonian momentum shell inte- 
gration. The form of ASk guarantees that the whole RG 
trajectory remains in the hypersurface of supersymmet- 
ric action functionals even within approximative calcula- 
tions; see 0, [l3| for further supersymmetric RG studies. 

In general, the full quantum effective action consists 
of all possible operators which are compatible with the 
symmetries of the theory. In this paper, we use a super- 
covariant derivative expansion to systematically classify 
these operators. Truncating the expansion at a given or- 
der constitutes a consistent approximation scheme. In a 
variety of bosonic and fermionic systems, such expansions 
have been proven to capture the physics of phase transi- 
tion and critical phenomena quantitatively, see, e.g. 



For our quantitative results presented below, we observe 
a satisfactory convergence: leading-order results generi- 
cally receive 0(1 — 10%) corrections at next-to- leading 
order, increasing up to ~ 50% for observables in the 
infinite-coupling limit. 



Let us start with the lowest-order derivative expansion, 
corresponding to a local-potential approximation for the 
superpotential. In this truncation, the effective action is 
approximated by 



Ffc = / dT 



(4) 



with a fc-dependcnt superpotential Wk- Projecting 
the flow equation ^ onto this truncation leads to 
the flow equation for the superpotential. For a sim- 
plified regulator choice with ri = and r2 = 
(|fc/p| — 1) 6* (l — (p^/A:^)), the superpotential flow re- 
duces to 



dtWk{<l>) = - 



ATI fc2 + FF^'(0)2 



(5) 



This regulator choice is convenient for analytical studies. 
We have verified explicitly by using a variety of differ- 
ent regulators that all conclusions drawn below from this 
simple choice also hold for other regulators. 

The corresponding flow equation for the dimensionless 
potential Wk — Wk/k admits a variety of nontrivial flxed- 
point solutions: already a polynomial expansion for Z2 
antisymmetric superpotentials to order n. 



dwk 



reveal n independent solutions to the fixed-point equa- 
tion dtWk = 0, in addition to the Gaufiian fixed point 20|. 
The latter has infinitely many relevant directions in 
agreement with perturbative power-counting in two di- 
mensions. The number of relevant RG directions de- 
creases for the other fixed points, culminating in a maxi- 
mally IR-stable fixed point with only one relevant direc- 
tion. This relevant direction exactly corresponds to the 
parameter in the superpotential, and it is common 
to all other fixed points as well. As is related to the 
bosonic mass term, we introduce a critical exponent vw 
for this direction, such that scales as ~ fc"^/""' at 
a fixed point. This critical exponent thus corresponds to 
the negative inverse of the associated eigenvalue of the 
stability matrix. The exponent uw plays a similar role 
for the superpotential, as the exponent v for the effective 
potential in Ising-type systems. At leading order in the 
derivative expansion, we obtain from Eq. ([5]): i^vk|lo = 1- 
Fixed-point potentials and critical exponents at next-to- 
leading order will be discussed below. 

In order to clarify the role of this relevant direction 
and its relation to supersymmetry breaking, we study 
the phase diagram of a Gaufiian Wess-Zumino model at 
infinite volume. This is defined in terms of a quadratic 
perturbation of the Gaufiian fixed point at the UV cut- 
off A: W'j^ — Aa(0^ — a\). The dimensionless control 
parameter 5 and coupling 7 can be associated with 



K^^Xko-I, 7:=2A-1Aa, 



(7) 
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FIG. 1: Top: Typical flow of the superpotential in the 
supersymmetry broken phase. Bottom: Phase diagram in 
the coupling-control-parameter plane for the Gaufiian Wess- 
Zumino model. The lowest set of data points corresponds to 
the lowest-order result ®. Higher orders converge rapidly. 



of higher orders of Eq. ([6]), we observe a rapid con- 
vergence of the phase boundary; see Fig. [1] (bottom). 
There exists a critical value for characterizing the 
phase transition in the strong coupling limit j —f oo. 
Our best estimate from a numerical higher-order solu- 
tion is Scr — 0.263; (the lowest-order result from Eq. ^ 
is ^7r/96 ~ 0.181). We conclude that supersymmetry 
can never be broken dynamically above this critical value. 
This agrees qualitatively with earlier results in the litera- 
ture 12 , 13, Q, [3] ■ Of course, the numerical value for 5cr 
is not universal but regulator-scheme dependent, simply 
because the quadratic initial conditions near the GauBian 
fixed point are nonuniversal. A quantitative comparison 
with other estimates of this critical value, e.g., taken from 
the lattice, is thus only meaningful if this scheme depen- 
dence is accounted for. 

The mass of the lowest bosonic excitation is given by 
the curvature of the full renormalized effective potential 
at the minimum. In the supersymmetric phase, fermionic 
and bosonic masses are identical and nonzero. In the bro- 
ken phase, the lowest fermionic excitation is a Goldstino 
with vanishing mass. The latter holds also at finite k 
in the supersymmetry-broken regime. As the flow in the 
IR is attracted by the maximally IR-stable fixed point, 
the bosonic mass in the fixed-point regime is governed by 
the only relevant direction a^. As the minimum of the 
effective potential in the broken regime is necessarily at 
a vanishing field, the bosonic mass yields 



where the scale is set by the UV cutoff A. In the follow- 
ing, we determine the region for the parameters (5, 7 for 
which supersymmetry is dynamically broken in the IR. 
A criterion for supersymmetry breaking is provided by 
a nonvanishing ground state energy, given by the min- 
imal value of V{(f>) = 5(W^fe^o(0))^- As a true phase 
transition only occurs in the infinite- volume limit [l^ . 
the ground state energy can only be identified at k —^ 0. 
Four snapshots of a typical flow of the superpotential in 
the broken phase are depicted in Fig. [1] (top) for S — 0.03 
and 7 = 0.2. 

For a simple estimate of the phase diagram, we trun- 
cate the superpotential at order 0^ in Eq. ([6|). From 
Eq. (O, we obtain the flow of the dimcnsionless couplings 



2fc^A 



k' 



(10) 
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2"^ 



6A"a 



2„2 



C\3 

a,A = -A+— , (8) 



which can be solved analytically. The phase transition 
curve in the (7, 5) plane is given by the initial values for 
which a^\k—,o changes sign. This condition yields 



arcsin(Q!) 



1 



TT 

3^- 



(9) 



Beyond this simple estimate, the phase transition curve 
can, of course, be identified by a full numerical inte- 
gration of the flow equation. Already upon inclusion 



where the latter proportionality holds in the vicinity of 
the IR fixed point where A — >const. and is governed 
by the critical exponent vyy- For vw > 1/2, the bosonic 
mass scales to zero upon attraction by the maximally IR- 
stable fixed point. We conclude that the broken phase 
remains massless in both bosonic and fermionic degrees 
of freedom. 

On the other hand, the limit fc ^ is an idealized IR 
limit. Any real experiment as well as any lattice simu- 
lation will involve an IR cutoff scale fcm characterizing 
the measurement, e.g., the scale of momentum trans- 
fer, detector size or lattice volume. Any measurement 
is therefore not sensitive to fc ^ but to fc — > fcm > 0. 
We conclude that a measurement of the bosonic mass in 
the broken phase will give a nonzero answer proportional 
to the measurement scale, whereas the goldstino will be 
truly massless. The bosonic mass in both phases as a 
function of the control parameter S is depicted in Fig. [2] 
(top) for a couphng of 7 = 0.2. 

At next-to-leading order in the derivative expansion, 
the fields acquire a wave function renormalization Zk, 
which is a multiplicative factor of the kinetic term in 
Eq. H]), £kin ^ Z|£kin- The flow of Zk gives rise to an 
anomalous dimension, 77 = —dt In Z^. Convergence of the 
derivative expansion requires rj < 0(1). 
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# nodes 


V 







0.4386 


1.2809 


1 


0.20 


1.11 


2 


0.12 


1.06 



TABLE I: Critical exponents of the first fixed points. 



From the polynomially expanded superpotential flow 
at next-to-leading order we can derive the following flow 
of the renormalized parameter a^: 



dtal 



2tt\ 4 



1 



'n\ 2 

2r* 



(11) 



At any fixed point with dt\ — > 0, we can read off that 
a?' again denotes a relevant direction with scaling a? ~ 
l^-i/fw where 



(12) 



This remarkable superscaling relation connects the su- 
perpotential exponent with the anomalous dimension 
J]. Recall that in Ising-like systems the thermodynamic 
main exponents (i.e., a, /3, 7 and d) are related among 
each other by scaling relations, and can be deduced from 
the correlation exponents u and r] by hyperscaling rela- 
tions. Beyond that, there is no general relation between 
ly and rj. The superscaling relation (|12p thus represents a 
special feature of the present supersymmetric model and 
is an exact relation to next-to-leading order in the super- 
covariant derivative expansion. Beyond next-to-leading 
order, Eq. p2|) might receive corrections, since higher- 
derivative operators involving the auxiliary field can still 
take influence on the flow of the superpotential. Never- 
theless, even if Eq. I|12p receives higher-order corrections, 
the most important result that a potentially more in- 
volved scaling relation among and rj exists still holds. 

Finally, we determine the possible nontrivial fixed- 
point superpotentials. Solving the coupled fixed-point 
equations for the dimensionless superpotential as well as 
the anomalous dimension self-consistently, we find a dis- 
crete set of fixed-point superpotentials with an increas- 
ing number of nodes, see Fig. [2] (bottom). The number 
of RG relevant directions at the fixed points in addition 
to the direction is equal to the number of nodes. The 
superpotential with no nodes is the next-to-leading or- 
der analogue of the maximally IR-stable fixed point dis- 
cussed above. Our quantitative next-to-leading order es- 
timates of the critical exponent vw and the anomalous 
dimension 77 for the fixed-point superpotentials of Fig. [2] 
(bottom) are summarized in Tab. |T1 As these quantities 
are universal, these results should constitute regulator- 
independent predictions. Of course, the truncation of the 
effective action introduces an artificial regulator depen- 
dence, the size of which can be used for an error estimate 
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FIG. 2: Top: Scale dependence of the renormalized bosonic 
mass. Bottom: Fixed-point superpotentials at next-to- 
leading order derivative expansion, x ~ Zk4>, u — w'l/Z]:. 



of our numerical results. From our regulator studies, we 
found a maximum variation of 10% on the positive criti- 
cal exponents. Subleading negative exponents can show 
larger variations and thus are difficult to estimate reli- 
ably, whereas the leading exponents vary even less, if at 
all [19]. 

Each fixed point constitutes a universality class and 
defines a distinct asymptotically safe UV completion of 
the Wess-Zumino model, i.e., RG trajectories emanating 
from different fixed points correspond to different theo- 
ries. Our findings are reminiscent of those in 2D scalar 
field theories 



17| , where the fixed points of the effec- 
tive potential can be related to conformal field theories 
[l8| . The connection between the present supersymmet- 
ric models at their fixed points and conformal field theo- 
ries remains an interesting question. 

To summarize, we have constructed a manifestly su- 
persymmetric flow equation for 2D M = 1 Wess-Zumino 
models. These models have nontrivial fixed-point super- 
potentials which can be classified by their relevant direc- 
tions. All fixed points share a relevant direction in the 
form of a bosonic mass term. At next-to-leading order, 
the associated critical exponent vw and the anomalous 
dimension r/ satisfy a superscaling relation which has no 
analogue in Ising-type spin systems. 

The fixed point of this relevant direction corresponds 
to the critical point separating the supersymmetric from 
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the broken phase. The initial condition for this rele- 
vant direction defines a control parameter, such that su- 
persymmetry breaking can be understood as a quantum 
phase transition. For the Gaufiian Wess-Zumino model, 
we observe that the control parameter stays finite even at 
arbitrarily large coupling in accord with a general argu- 
ment by Witten. In the broken phase, our superscaling 
relation predicts that the measured bosonic mass is pro- 
portional to the momentum scale set by the detector. 

Helpful discussions with G. Bergner, C. Wozar, T. Fis- 
chbacher and T. Kaestner are gratefully acknowledged. 
This work has been supported by the Studienstiftung 
des deutschen Volkes, and the DFG under GRK 1523, 
Wi 777/10-1, FOR 723 and Gi 328/5-1. 



[1] E. Witten, Nucl. Phys. B188, 513 (1981). 

[2] E. Zhao and A. Paramekanti, Phys. Rev. Lett. 97, 230404 

(2006) . 

[3] B. Uchoa and A. H. C. Neto, Phys. Rev. Lett. 98, 146801 

(2007) . 

[4] J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 



[5] A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

[6] P. Strack, S. Takei, and W. Metzner (2009), 

arXiv:0905.3894 [cond-mat.str-el]. 
[7] C. Wetterich, Phys. Lett. B301, 90 (1993). 
[8] F. Synatschke, G. Bergner, H. Gies, and A. Wipf, JHEP 

03, 028 (2009). 
[9] H. Sonoda and K. Ulker (2008), arXiv:0804.1072. 
[10] O. J. Rosten (2008), arXiv:0808.2150. 
[11] B. Delamotte (2007), arXiv:cond-mat/0702365. 
[12] E. Witten, Nucl. Phys. B202, 253 (1982). 
[13] M. Beccaria, G. F. De Angelis, M. Campostrini, and 

A. Feo, Phys. Rev. D70, 035011 (2004). 
[14] M. Beccaria, M. Campostrini, and A. Feo, Phys. Rev. 

D69, 095010 (2004). 
[15] J. Ranft and A. Schiller, Phys. Lett. B 138, 166 (1984). 
[16] T. R. Morris, Phys. Lett. B345, 139 (1995). 
[17] R. Neves, Y. Kubyshin, and R. Potting (1998), hep- 
th/9811151. 

[18] A. Zamolodchikov, Yad. Fiz. 44, 529 (1986). 

[19] F. Synatschke, H. Gies, and A. Wipf (2009), 
arXiv:0907.4229 [hep-th]. 

[20] Beyond the polynomial expansion, we find fixed-point so- 
lutions corresponding to oscillating sine-Gordon-type po- 
tentials, bouncing potentials as well as potentials which 
confine the field to a compact target space; see [l^ for a 
detailed discussion. 



